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Abstract 



The exact solution of the asymmetric exclusion problem with N distinct 

classes of particles (c = 1, 2, . . . , A'^), with hierarchical order is presented. In 

this model the particles (size 1) are located at lattice points, and diffuse with 

equal asymmetric rates, but particles in a class c do not distinguish those in 

O ■ the classes c > c from holes (empty sites). We generalize and solve exactly 

O 

this model by considering the molecules in each distinct class c = 1,2, . . . ,N 
j^ ' with sizes Sc {sc = 0,1,2,...), in units of lattice spacing. The solution is 

derived by a Bethe ansatz of nested type. 
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I. INTRODUCTION 

The similarity between the master equation describing time fluctuations in nonequihb- 
rium problems and the Schrodinger equation describing the quantum fluctuations of quan- 
tum spin chains turns out to be fruitful for both areas of research |T|]- ||T5|. Since many 
quantum chains are known to be exactly integrable through the Bethe ansatz, this provides 
exact information on the related stochastic model. At the same time the classical physical 
intuition and probabilistic methods successfully applied to nonequilibrium systems give new 
insights into the understanding of the physical and algebraic properties of quantum chains. 

An example of this fruitful interchange is the problem of asymmetric diffusion of hard- 
core particles on the one dimensional lattice ( see [p!^-[18|| for reviews). This model is related 



to the exactly integrable anisotropic Heisenberg chain in its ferromagnetic regime [|T^ (XXZ 
model). However if we demand this quantum chain to be invariant under a quantum group 
symmetry Ug{SU{2)), we have to introduce, for the equilibrium statistical system, unusual 
surface terms, which on the other hand have a nice and simple interpretation for the related 
stochastic system ^j^. 

In the area of exactly integrable models it is well known that one of the possible extensions 
of the spin-i XXZ chain to higher spins, is the anisotropic spin-S Sutherland model (grading 
ei = €2 = ■ ■ ■ = e2s+i = 1) [^. On the other hand in the area of diffusion limited reactions 
a simple extension of the asymmetric diffusion problem is the problem of diffusion with 



particles belonging to N distinct classes (c = 1,2,..., N) with hierarchical order |22|- ||2^ . 
In this problem a mixture of hard-core particles diffuses on the lattice. Particles belonging 
to a class c (c = 1, . . . , A^) ignore the presence of those in classes c' > c, i. e., they see 



them in the same way as they see the holes (empty sites). In it was shown that for open 
boundary conditions the anisotropic spin-1 Sutherland model and this last stochastic model, 
in the case N = 2, are exactly related. The Hamiltonian governing the quantum or time 
fluctuations of both models being given in terms of generators of a Hecke algebra, invariant 
under the quantum group UgSU{3). In fact this relation can be extended to arbitrary values 
of A^, and the quantum chain associated to the stochastic model is invariant under the 
quantum Ug{SU{N + 1)) group. In this paper we are going to derive through the Bethe 
ansatz the exact solution of the associated quantum chain, in a closed lattice. Recently |1^ 



(see also |T^) we have shown that without losing its exact integrability, we can consider 
the problem of asymmetric diffusion with an arbitrary mixture of molecules with different 
sizes (even zero) , as long they do not interchange positions, that is, there is no reactions. In 
this paper we are going to extend the asymmetric diffusion problem with A^ type of particle 
with hierarchical order, to the case where the particles in each class have an arbitrary size, 
in units of the lattice spacing. Unlike the case of asymmetric diffusion problem, we have in 
this case a nested Bethe ansatz [^]. A pedagogical presentation for the simplest case N = 2 



was presented in [p6 |. 



The paper is organized as follows. In the next section we introduce the generalized 
asymmetric model with A^ type of particle with hierarchical order and derive the associated 
quantum chain. In section 3 the Bethe ansatz solution of the model is presented. Finally 
in section 4 we present our conclusions, with some possible generalizations of the stochastic 
problem considered in this paper, and some perspectives on future work. 



II. THE GENERALIZED ASYMMETRIC DIFFUSION MODEL WITH N 
CLASSES OF PARTICLES WITH HIERARCHICAL ORDER 

A simple extension of the asymmetric exclusion model, where hard-core particles diffuse 
on the lattice, is the problem where a mixture of particles belonging to different classes 
(c = 1, 2, . . . , N) diffuses on the lattice. This problem in the case where we have only N = 2 
classes was used to describe shocks |^2|- |^^ in nonequilibrium and also has a stationary 



probability distribution that can be expressed via the matrix-product ansatz P7[. In EH 



it was also shown that the stationary state of the case A^ = 3 can also be expressed by 
the matrix-product ansatz. In this model we have ni,n2, ■ ■ -Un molecules belonging to the 
classes c = 1, 2, . . . , iV, respectively. All classes of molecules diffuse asymmetrically, but 
with the same asymmetrical rates, whenever they encounter empty sites (holes) at nearest- 
neighbor sites. However, when molecules of different classes, c and c' (c < c'), are at their 
minimum separation, the molecules of class c exchange position with the same rate as they 
diffuse, and consequently the molecules in the class c see no difference between molecules 
belonging to the classes c' > c and holes. 

We now introduce a generalization of the above model, where instead of having unit size, 
the molecules in each distinct class c = 1, 2, . . . , A^ have in general distinct sizes si, S2, . . . , sat 
(si, . . . , stv = 1, 2, . . .), respectively, in units of lattice spacing. In Fig. 1 we show some 
examples of molecules of different sizes. We may think of a molecule of size s as formed by 
s monomers (size 1), and for simplicity, we define the position of the molecule as the center 
of its leftmost monomer. The molecules have a hard-core repulsion: the minimum distance 
(io/3, in units of the lattice spacing, between molecules a and /3, with a in the left, is given 



by dai3 = Sa- In order to describe the occupancy of a given configuration of molecules we 
attach at each lattice site i {i = 1,2, . . . ,L) a variable /3j {i = 1,2, ... ,L), taking the values 
Pi = 0,1, ... ,N. The values /3 = 1,2, ... ,N represent the sites occupied by molecules of 
class c = 1,2, ... ,N, respectively. On the other hand the value f3 = represent the empty 
sites as well the excluded ones, due to the finite size of the molecules. As an example, in 
a L = 8 sites chain, the configuration where a particle of class 1, with size si = 2 located 
at site 1, and another particle of class 2, with size S2 = 3 located at site 3, is represented 
by {P} = {1,0,0,2,0,0,0,0}. Then the allowed configurations are given by the set {/?j} 
(i = 1, . . . , L), where for each pair {Pi, Pj) 7^ with j > i we should have j — i > sp.. 

The time evolution of the probability distribution P{{P}, t), of a given configuration {/?} 
is given by the master equation 

^^^ = E i-nm - m)Pim, t) + mp'} -. {p})p{{p'}, t), ] (i) 

aqui where r({/9} -^ {P'}) is the transition rate for configuration {P} to change to {P'}. In 
the present model we only allow, whenever the constraint of excluded volume is satisfied, the 
particles to diffuse to nearest-neighbor sites, or to exchange positions. The possible motions 
are diffusion to the right 

A 0.+1 ^ 0. A+i, (/3 = 1, . . . , iV) (rate Tr) (2) 

diffusion to the left 

0, A+i ^ P^ 0^+1, iP=l,...,N) (rate T^) (3) 

and interchange of particles 

A P:+s, ^ A' P^+s,, , {P<P' = 1,...,N) (rate Tr) 



A A'+s, ^ A' A+.,„ (/3 > A = 1, • • • , iV) (rate T^). (4) 

As we see from (^ , particles belonging to a given class c interchange positions with those of 
class d > c with the same rate as they interchange positions with the empty sites (diffusion). 
We should remark however that unless the particles in class c' have unit size {sc = 1), the 
net effect of these particles in those of class c is distinct from the effect produced by the 
holes, since as the result of the exchange the particles in class c will move by s^' lattice size 
units, accelerating its diffusion. 

The master equation (||) can be written as a Schrodinger equation in Euclidean time (see 
Ref. for general apphcation for two body processes) 

if we interpret \P >= P({/3}, t) as the associated wave function. If we represent A as |/3 >i 
the vector |/? >i ®|/9 >2 ® ■ ■ ■ ® |/3 >l will give us the associated Hilbert space. The process 
(I)- (I) gives us the Hamiltonian (see Ref. for general applications) 



j 

N 

a=l 

N N 



with 



D = Tn + TL, 6+ = -!^, e^ = —^ (^^ + ^^ = 1), (7) 

' e+ a < /3 

ea/3 = <^ a = /3 (8) 

e_ a > P 
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and periodic boundary conditions. The matrices E""'^ are (A^ + 1) x (A^ + 1) matrices with 
a single nonzero element {E°''^)ij = 5a^i5pj {a,P,i,j = 0, . . . , A^). The projector V in (^, 
projects out from the associated Hilbert space the vectors \{P} > which represent forbidden 
positions of the molecules due to their finite size, which mathematically means that for all 
i,j with Pi, Pj 7^ 0, \i — j\ > S/^^ (j > i). The constant D in (^ fixes the time scale and 
for simplicity we chose D = 1. A particular simplification of (P) occurs when the molecules 
in all classes have the same size Si = S2 = ■ ■ ■ = sn = s. In this case the Hamiltonian 
can be expressed as an anisotropic nearest-neighbor interaction spin-A^/2 SU{N + 1) chain. 
Moreover in the case where their sizes are unity [s = 1) the model can be related to the 



anisotropic version [^ of the SU{N + 1) Sutherland model |Q with twisted boundary 
conditions. 

III. THE BETHE ANSATZ EQUATIONS 
We present in this section the exact solution of the general quantum chain (0). A 



pedagogical presentation for the particular case where N = 2 was presented in [^ . 

Due to the conservation of particles in the diffusion and interchange processes the total 
number of particles rii, n2, . . . , n^ in each class are good quantum numbers and consequently 
we can split the associated Hilbert space into block disjoint sectors labeled by the numbers 
rii, n2 ■ ■ ■ riN (n^ = 0, 1, . . . ; i = 1, . . . , iV). We therefore consider the eigenvalue equation 

if |ni, n2, . . . , riAT >= E\ni, 77-2 . . . , wat >, (9) 

where 

|ni, n2, . . . , riAT >= ^ ^ /(xi, Qi; . . . ; x„, Qn)\xi, Qi; . . . ; x„, Q„ >, (10) 

{Q} {^} 



and n = J2iLi ^i is the total number of particles. In (|ToD |xi, Qi; . . . ; a;„, Q„ > means the 
configuration where a particle of class Qi {Qi = 1,2, ... , N) is at position Xi {xi = 1, . . . ,L). 
The summation {Q} = {Qi, . . . , Qn} extends over all permutations of the n integer numbers 
{1, 2, ... , N} in which n, terms have the value i {i = 1, 2, ... , N) , while the summation 
{x} = {xi, . . . , Xn} runs, for each permutation {Q}, in the set of the n nondecreasing integers 
satisfying 

x^+l>Xi + SQ^, i = l,...,n-l, Sq^< Xn-Xi< N - Sq„. (11) 

Before getting the results for general values of n let us consider initially the cases where we 
have 1 or 2 particles. 

n = 1. For one particle on the chain, in any class c = 1, 2, . . . , A^, as a consequence of 
the translational invariance of (§) it is simple to verify directly that the eigenfunctions are 
the momentum-fc eigenfunctions 

L 

\0,...,0,l„0,...,0>=J2f{x,c)\x,c>, c = l,...,N (12) 



x=l 

with 

Ovr/ 
f{x,c)=e''\ k = —, / = 0,1,...,L-1, (13) 

and energy given by 

E = e{k) = -{t_e''' + t+e-''' -I). (14) 

n =2. For two particles of classes Qi and Q2 {Qi, (^2 = 1, 2, . . . , A^) on the lattice, the 
eigenvalue equation (|^) gives us two distinct relations depending on the relative location of 
the particles. The first relation applies to the case in which a particle of class Qi (size sqJ 

8 



is at position xi and a particle Q2 (size sq^) is at position X2, where X2 > xi + sq-^. We 
obtain in this case the relation 

Ef{xi, Qi; X2, Q2) = -e+f{xi - 1, Qi; Xs, Q2) - ^-f{xi, Qi; X2 + 1, Q2) 
-e-f{xi + 1, Qi] X2, Q2) - e+/(a;i, Qi]X2- 1, (^2) + 2/(xi, Qi; X2, (^2), (15) 

where we have used the relation e+ + e_ = 1. This last equation can be solved promptly by 
the ansatz 

/(xi,Qi;x2,Q2) = ^A^^^;g^e*(^'^i^^+^'^2^^) 
p 

__ aQi,Q2 i{k\xi+k2X2) _|_ /\Qi,Q2 i{k2Xi+kiX2) (IS) 

with energy 

E = e{h) + e{k2), (17) 

where ki,k2,Ai2 and ^21 are free parameters to be fixed. In (|1^) the summation is 



over the permutations P = Pi, P2 of (1,2). The second relation applies when 0:2 = xi + sqj. 



In this case instead of ([T5| ) we have 



Ef{xi, Qi]Xi + sq^,Q2) = -e+/(xi - 1, Qi, Xi + sq^, Q2) - e_/(xi, Qi; Xi + sq^ + 1, Q2) 
-^Q2,Qif{xi, Q2; xi + sq„ Qi) + (1 + eQi,Q2)/(a;i, Qi, Xi + sg,, ^2)- (18) 

If we now substitute the ansatz ( P^B] ) with the energy (jT^), the constants A12 and A21 , 
initially arbitrary, should now satisfy 

E{ f^Pi,P2 + e^'^^ (1 - 6Q„Q,)1 e^'^-2(^Q.-i)A«-«/ + iQ2,Q.e''^^'^^A^p^;g^} = (19) 



p 



where 



n 



Lrri 



(e+ + e_e' 



i(ki+km)-\ 



(20) 



At this point it is convenient to consider separately the case where Qi = Q2 from those 
where Q, ^ Q^. li Q, = Q^ = Q (Q = 1, . . . , N) eq. (0) gives 



p 



(Vp.^P, + e'^^^) e'^P2i.'Q-^) Af^% = 



(21) 



and the cases Qi 7^ Q2 give us the equations 



E 
p 



Akp 



gifcp2(sQ^-l)^Ql,Q2 



gJfcp2(sQ2-l)^Q2,Ql 



0. 



Performing the above summation we obtain, after lengthy but straightforward algebra, the 
following relation among the amplitudes 



^Q2,QlgJfc2(sQ2-l) 



ik\ 



X < 



Pi,2 + e 

" Pi 2 + 6*^^2 



l-$(fcl,fc2) 






^Qi,Q2gifci(sQj-l) 

^Q2,Qlgifcl(sQ2-l) 



where 



$(fci, A;2) 



gifci _ gifc2 



(22) 



Equations ( pTD and (^4 ) can be written in a compact form 



^Pl,P2 



TV 



-Pl,P2 Z^ 



5g/;g,^(A;p„ A;pJA?^;gS (Qi, g2 = 1, • • • , iV) 



(23) 



with 



-«j 






^tkj 



(24) 
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where we have introduced the 5* matrix. From ( pi]) and (p2|) this S matrix has only A^(2A^— 1) 
non zero elements, namely 

SSlSliku k2) = [1 - eQ.^QMh, h)] e^('=-^^)(^Qi-i) (Qi, Q^ = 1, • • • , iV), 

SqlfA^i. h) = eQ,,QMku fc2)e^'^(^«^-^)e-^'=^(^«i-i) {Q,,Q, = 1, . . . , N;Q, ^ Q2). (25) 

Equations (^) do not fix the "wave numbers" ki and k2. In general, these numbers are 
complex, and are fixed due to the cychc boundary condition 

fixi, Qi; X2, Q2) = fix2, Q2; xi + N, Qi), (26) 

which from (0) give the relations 

AQ1Q2 _ ^ikxN /\Q2,Qi aQi,Q2 _ ^ik2N aQ2,Qi (07\ 

/±12 —^ ^2,1 ) ^2,1 —^ ^2,1 • \^') 



This last equation, when solved by exploiting (|23D -(P5D, gives us the possible values of ki 
and ^2, and from ([TTD the eigenenergies in the sector with 2 particles. Instead of solving 
these equations for the particular case n = 2 let us now consider the case of general n. 

General n. The above calculation can be generalized for arbitrary occupation 
{rii, n2, . . . , TiAf} of particles in classes 1, 2, . . . , A^, respectively. The ansatz for the wave 
function ( [T0| ) becomes 

/(xi, Qi; . . . ; x„, Q„) = ^ ^g- J;e^('=-i-i+-+'=-"^'"), (28) 

p 

where the sum extends over all permutations P of the integers 1, 2, . . . , n, and n = J2iLi ^i 

is the total number of particles. 

The application of the translation operator in the above wave functions implies that (|TUp 

are also eigenf unctions of the momentum operator with eigenvalues 

11 



" 27r/ 

P = T.k, = ^, (/ = 0,1,...,L-1). (29) 

For the components \xi, Qi, . . . ; x„, Qn > where Xj+i — Xj > sq- for i = 1,2, ... ,n, it is 
simple to see that the eigenvalue equation (|) is satisfied by the ansatz (pS]) with energy 

n 

E = Y.e{k,). (30) 

On the other hand if a pair of particles of class Qi, Qi+i is at positions Xi, Xj+i, where Xj+i = 
Xj + sq., equation (|) with the ansatz (|^) and the relation (|30D give us the generalization 
of relation (^31), namely 



N , , 



;iiV2 



(31) 



with 5* given by eq. (25). Inserting the ansatz ( psD in the boundary condition 

/(xi, Qi;...;xn, Qn) = f{x2, Q2]---] Xn, Qn] Xl + N, Qi) (32) 

we obtain the additional relation 

AQl,--,Qn ikp.N AQ2,---,Qn,Ql /"qq^ 

^Pl,...,Pn - ^ ^P2,...,P„,Pl ' W"^J 



which together with (pT]) should give us the energies. 



Successive applications of (|3TD give us in general distinct relations between the am- 
plitudes. For example ^''''^^ ^j'^'g,... relate to ^^^^'^3 ^j'^j,... by performing the permutations 

a/37 ~^ P^^l ~^ Pl^^ ~^ iP^^ o^^ (^Pl ~^ (^iP ~^ l^^P ~^ iPciy and consequently the S'-matrix 



should satisfy the Yang-Baxter [T^J2^ equation 

AT 

J2 S"',Y ih , k2)SJly',{ki , k3)Sj',']i!,{k2 , k^) = 

E s';>:;:{k2,k,)s"j:{k,,k,)s]^];,{k,,k2), (34) 

7,7',7"=1 

12 



for a,a',a", P, P', P" = 1,2,...,N and 5* given by (]25|). Actually the relation (0) is a 
necessary and sufficient condition [0,^ to obtain a non-trivial solution for the amplitudes 



in Eq. {^. 

We can verify by a long and straightforward calculation that for arbitrary number of 
classes N and values of the sizes Si, S2, ■ ■ ■ , s^, the S matrix (P^), satisfies the Yang-Baxter 
equation (|5), and consequently we may use relations (^TJ) and (^) to obtain the eigenen- 
ergies of the Hamiltonian (|^). Applying relation (|3lD n times on the right of equation (|33D 
we obtain a relation between the amplitudes with the same ordering in the lower indices: 

AQl,...,Qn _ Akp^N AQ2,--;Qn„Ql _ ( TT ■^ \ Jkp^ N \^ \^ 

V»=2 / q;,...,qj,q'/,...,q^ 

•^Q'l ,'q'/ ^ ^^1 ' ^^1 ^ ^Q'2 .'qI \^p-2^^Pi)' ' ' ^qI_ ^ \q'1_ ^ ( ^Pn- 1 5 ^Pi ) 'S'q; ,'qA ( kp^ , /cp^ ) Apl'^ pj' , ( 35 ) 

where we have introduced the harmless extra sum 

N ^ „ 

1= Y. ^Q'lQ[^Q'{,Qi = J2 Sgl^^^^ikp^kpJ (36) 



see 



26[] for ilustrations of the above equations). In order to fix the values of {kj} we should 



solve (|35D, i.e., we should find the eigenvalues A{k) of the matrix 



nk)\^\= E (l[SQ^^f{kp^,k)], (37) 

Q'/,...,Q^=l \/=l / 

with periodic boundary condition 

The Bethe-ansatz equations which fix the set {ki} will be given from (|35D by 

e-''^'' = {-ir'{U.^i,j^Hkj), J = l,...,n. (39) 

13 



The matrix T{k) has dimension A^" x A^" and can be interpreted as the transfer matrix of a 
inhomogeneous A^(2A^ — l)-vertex model in a two dimensional lattice with periodic boundary 
condition in the horizontal direction (n sites). Due to the special form of the S matrix (pSj) 
the eigenvalues of ( |57D are invariant under a local gauge transformation where for each factor 
S{kp^,k) in i^: 

S^^^{kp,k) -. S^^^{kp,k)^, (40) 

where 0^'^ (/ = 1, . . . , L; a = 1, . . . , A'") are arbitrary functions. If we perform the trans- 



formation (^OD with the special choice 



e 



ikp^(s^~i)^ (/ = 2,3,...,A^), (41) 



0« 
the equivalent transfer matrix to be diagonalized is given by 

f (fc){« = e-'=Sr..(^«»-)To(fc)|« (42) 

where 

N ( "^ ^n nil \ 

ro(fc)g.\= E n^|Qf(A:p,,A:) , (43) 

Q'i...,Q':=\ \i=\ / 

with the twisted boundary condition 

St3t\kp.. k) = S^r^{kp„, k)<l>Qii (44) 



with twisted phase 



$^ = e^(-.-i)E;=i'=.^ / = !,..., AT. (45) 



The matrix S in (^) and (^) is obtained from those in ( pS] ) by taking the size of all 
particles equal to the unity. In this way the problem is transformed into the evaluation of 

14 



the eigenvalues of a regular (all particles with size 1) inhomogeneous transfer matrix Tq with 
n{2N — 1) non-zero vertex and twisted boundary condition. 
Diagonalization of To(/c) 

The simplest way to diagonalize Tq is through the introduction of the monodromy matrix 
A4{k) |]25|, which is a transfer matrix of the inhomogeneous vertex model under considera- 
tion, where the first and last link in the horizontal direction are fixed to the values /ii and 
/in.+i (/Ui, /i„+i = 1, 2, . . . , A^), that is 

N 

The monodromy matrix -M \qi\^!^"^^ {k) has coordinates {Q}, {Q'} in the vertical space (A^" 
dimensions) and coordinates fii, fin+i in the horizontal space (A^^ dimensions). This matrix 
satisfy the following important relations 

N 

N 

E -^!:},:r^(^') -^l';}::r^(^)^:;:x^(^''^)' (47) 

for fii, z/i, /i„+i, Un+i = 1, 2, . . . , A^.This relation follows directly from successive applications 



of the Yang-Baxter equations ( ^4]) (see ||26|, for a graphical representation of these equations). 
In order to exploit relation ( ^7]) let us denote the components of the monodromy matrix 
in the horizontal space by 

c{kr = Mi::^^^jk), D{k) = Mi::}^^^{k), (48) 
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where (a,/3 = 1, 2, . . . , A^ — 1). Clearly the transfer matrix To{k) of the inhomogeneous 
lattice with twisted boundary conditions, we want to diagonalize, is given by 

N-l 

Uk)=Y.K{k) + D{k). (49) 

As a consequence of (47) the matrices A", 5", C" and D in (|48| ) obey some algebraic 
relations. By setting {vi, ^i.Vn+i^ l^n+i) = {.N,a,'~^, [3) in (^Tf ) we obtain 

Af(fc)i^^(fcO = - ^y.';> ^(A:)AZ(A:0+ E ^S0^^"'(^')< W, (50) 

^N,a['^j'^) a',l3'=l ^N,a[^ ^^) 

with (a, /3 = 1, . . . , A^ — 1). By setting (t/i, /xi, z/„+i, yUn+i) = {N, N, N, a) we obtain 

Dik)B-ik') = ^$/^B'^ik')Dik) - |sMls"(A:)D(fc'), (51) 



where (a = 1, . . . , A^ — 1). The diagonalization of To{k) in (^) will be done by exploiting 
the above relations. This procedure is known in the literature as the algebraic Bethe ansatz 
. The first step in this method follows from the identification of a reference state \^ >, 



which should be an eigenstate of A'j^{k) and D{k), and hence To{k), but not of B"'{k). In 
the present case a suitable reference state is \Q >= \{ai = N} >;=!,...,„, which corresponds 
to a state with N-class particles only. It is simple to calculate 

A^^{k)\Q> =a^^{k)\Q>, D{k)\Q>=d{k)\Q>, 

n 

C^ik)\n > = 0, S"(A;)|fi >= E&r(A:)|fi« >, (52) 



i=l 



where 



i=l i=l 

i— 1 n 

bm = $^ n S^lNikPnk) n SZikp^k), (53) 

1=1 l=i 
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and 1^2^ >= |{q;^_^j = N},ai = a >. The matrices B°'{k) act as creation operators in the 
reference ("vacuum") state, by creating particles of class a (1,2,..., N) in a sea of particles 
of Nth class \VL >. We then expect that the eigenvectors of Toi^k) corresponding to nii 
(1,2,..., n) particles, belonging to classes distinct from A^, can be expressed as 

m 

where {k\ , / = 1, • • • , rrii} and -F/3i,...,/3„ are variables to be fixed by the eigenvalue equation 

To{k)\kl^\ F >= M°\k)\kl^\F> . (55) 

Using (|50D successively, and (|5^) , (|53|) we obtain 

N N 

A:{k)B^^{ki'^)B^^{k?) . ■■B^-^{k^^l) = E E 

K,...,<H=l}{/5i.-./3mi=l} 
c/3i,a (uW l.\q^^'^'ui.W l.\ Q^'"l-l'^mi-2/, (1) ,X r.^"H'/3mi-i/, (1) ,-. 

x^^ '=\y\': ' B^ik'l^^B'^'^ik^^^) ■ • •5<i(A;(i|)|fi > + "unwanted terms", (56) 

\.lj=l'-'N,a\'^j 5 '^J 

where the "unwanted terms" are those ones which are not expressed in the "Bethe basis" 
produced by the B"{kj ) operators. Similarly, using (^TJ) successively and (^)-(|5B|) we 
obtain 



Dik)Bf'^ik\'^)...B^-^ik^^l)\n>=^r,[l[SN:Nihk) 



X iTT 4'^ m U^nfe'^)---^^'"Kfe^^)l^> + "unwanted terms". (57) 



The relations (|56D and (^) when used in (|54)-(|55D give us 



qN,l 



n"=i^^jv,'i(fcj,fc) y- y- T^ n.^{/3} F. „, R"'i ^1,(1)^ 

llj=l'^Ar,ll'Cj ,«;; {/3}{a'} 



mi cA^.V^Cl) 
mi 

+$^ n S^lNih. A:) E n I f^7rilT)T 1 FmB'^ik^^') " " " ^'^'"^ (^S)!^ > 
i=l {/3}«=i WftA^'^i )/ 

+" unwanted terms", (58) 
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where 

N-l N-1 /mi-1 \ fl fl' 

Tmfj, = E $. E 4is(^s^A:) n ^5i.(^J^^) ^.:;^'^(e>^) (59) 

a=i /3;,...,/3;,^=i \ i=i ' ' ) 

is a (A^ — 1)™^ -dimensional transfer matrix of a inhomogeneous vertex model, with inhomo- 
geneties {A;^|, femi-i5 • • • ; ^i } (notice the reverse order of the inohomogeneties, when com- 
pared with (PSP) and twisted boundary conditions (boundary phases $q,, a = 1, . . . , A^ — 1). 
In order to proceed we need now to diagonalize the new transfer matrix Ti[k), that is 
we must solve 

E^l(^){a'}^{/3}=A«(fc)^K} (60) 

m 

and then (|58D give us 

To{k)\k^^\ F >= A'^^\k)\k^^\ F > + "unwanted terms", (61) 

where, using the fact that 5'^'jy (/c^, /c) = 1, 



^'°'W-i!r|^SS^^"W + *"nwT7fTaTT- m 



I\j=i S^'^i [kj , k) 1=1 Sj^'^i [k, ki 

In order to proof that A*^*^^ and \K2'\ F > are the eigenvalues and eigenvectors of To(A;), we 
should fix {k[ , . . . , ^mj} by requiring that the "unwanted terms" in (^) vanish. Although 



for N = 2 this calculation is not complicated PB[ for arbitrary A^ it is not simple. Since the 



expression (|6^) for the eigenvalues should be valid for arbitrary values of k we can obtain 



A^^\kj ) in an alternative way from the following trick |3TI]. At fc = A;] (j = 1, . . . ,mi 



the denominators of the factors in ( |62D vanishes (S^J (A;]- ,kj) = 0,l^N), and since we 
should have a finite result, we have the conditions 



1 mi qN,U,{1) ,(1)x 



A")(fcf)^t.n ,.a„ ,.!.> n ss :;;,•,' . j=i.. ....... (63) 



Notice that our result in ( |63D does not depend on the particular ordering of the additional 
variables kj (j = 1, . . . , mi). This means that if instead of the ordering chosen in (0), we 
chose the reverse order, namely, 

\k^a\F>=Y.F,,,...,,^^B^-^iki:;l)B^-^-^ikill,) ■ ■ ■ B\k^^^)\n > (64) 

1/3} 



we obtain the same results (|6lD -(|63|) but now Ti is the transfer matrix, with boundary 
condition specified by the phase $„, of a problem with (iV — 1) species and inhomogeneities 
k\ , . . . , fc^l (notice we have now the same order of the inhomogeneties as in (^3]) ) . This 
means that the eigenvalue A(A;) = A^^\k) of the transfer matrix of the problem with A^ 
classes and inhomogeneities {k[ ,k2 , ■ ■ ■ , k^'>) = (fci, ^2, . . . , kn) is related to the eigenvalue 
A*^^)(A;) of the problem with (A^ — 1) classes and inhomogeneities /cj , /ca ? . . . , k':^ . Iterating 
these calculations we obtain the generalization of the relation (p^ ) and the condition (|63| ) 



/ rni \ mi+i , \ 

A'"(*) = n s?:;(4", *) n ..,,1., ,, a"«>(*) + 

\i'=i ) \v=\ S^\{k\, \k) 



^""^ ,.,.,,, .IV ^ = 0.1.- --.A^-l. (65) 

l'=l D^i l^K, Ki, ) 






A"^"(>=r')-<'^-.n L«.i„. ,mK n .^;1.;1.. . (66) 



which connects the eigenvalues of the inhomogeneous transfer matrix Ti{k) and T/+i(A;), with 
inhomogeneities {kj } and {kj }, related with the problem with (A^ — /) and (A^ — / — 1) 
classes of particles, respectively. 

However from ( |39D and (^)-(^), in order to obtain the Bethe-ansatz equations for our 
original problem we need the eigenvalues of the transfer matrices evaluated at kj {j = 
1, . . . , n), i. e., A^'^\kj), which are given by 
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mi T 

A^°HA:pJ = $.n ..,w, ,av (67) 

since 11^=1 S^'aikj, kp-) = 0. The conditions that fix the variables {kj , j = 1, . . . ,mi) are 
given by (|63|). In the left side of this equation we have A^^\kj ), which are the eigenvalues 
of the transfer matrix Ti of the model with (A^ — 1) classes of particles and inhomogeneities 
{kj \i = 1, . . . , mi}, evaluated at the partcular point k^ . This value can be obtained from 
(|55|) which gives a generalization of ( p7| ) 



m-i+i 



1 



A^')(fcf) = ^^-^n .^,vJo ,a+i), (^ = o,i,...,iv-i). (68) 

i'=l >Jiv^i(,ftj , ft;' J 



The condition (|63[) is then replaced by 

A'"(4") - ^A'-i n .«..,,., , ..I, = -^^ n ..,,,,., ,a,, n .":;;.)',;.) ■ («"' 

Z'=l Ojv^j^ l,ftj , ft/' } i=lD]^]^{K,^ , Kj ) l'=i^l'j^j Op^i^Kj , % ) 

where now we need to find the relations that fix {k, }. Iterating this process we find the 
generalization of (BO) 



mi+i -, mi_i 



1 






' '"J 



^ n e^i Jo'jo '(^' = ^'^---^'' / = 0,l,...,iV-2). (70) 

Equations (|67|) and (|70D give us the eigenvalues of the tranfer matrix To(/c) evaluated at the 
points {fcj}, i. e. A.^'^\kj). Inserting the above results in (^21) and then in (|39D we obtain 
the Bethe-ansatz equations of our original problem. 

The eigenenergies of the Hamiltonian (^ in the sector containing rii particles in class i 
(i = 1,2, . . . , N) [n = Z^jLi f^j) and total momentum p = ^ (/ = 0, 1, . . . , L — 1) are given 
by 



" (0) _ii.m 



Y^/, Jky' , , „-ik 



E = ~y:{e.e''. +e^e-''^ -1), (71) 
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where {kj ' = kj,j = 1, . . . , n} are obtained from the solutions {fc^ , Z = 0, . . . , A^ — 1; j 
1, . . . , m;} of the Bethe ansatz equations 






X 



- 6+(e<'-e^^r) 



n ^nT)-(5^ T(5) J = l,2,...,n, (72) 






and 



13=1 e+ + e_e^ " z' ' - e f s=i e+ + e.e*^"'? +'=« ^ - e*'=<^ 



X 



e+ + e_e^ " <^' ' - e 



n ,,.a+i),,a+i). T(ITT) / = 0,l,...,iV-2; a = l,...,m;, (73) 

and m; = Y.f=i f^j, ^ = 0,. . . ,N (mo = N^mjq = 0). It is interesting to observe that in 
the particular case where n2 = n^ = ■ ■ ■ = n^ = we obtain the Bethe-ansatz equations, 
recently derived |T5[ (see also [1^), for the asymmetric diffusion problem with particles of size 
si. Also the case si = S2 = ■ ■ ■ = sn = ^ give us the corresponding Bethe ansatz equations 
for the standard problem of A^ type of particles in hierarchical order. The Bethe-ansatz 
solution in the particular case of N=2 with a a single particle of class 2 (ni = n — l,n2 = I) 
was derived recently [^. The Bethe-ansatz equations for the fully asymmetric problem is 



obtained by setting in (|72D-(|73D e+ = 1 and e_ = 0. 

IV. CONCLUSIONS AND GENERALIZATIONS 

We obtained through the Bethe ansatz the exact solution of the problem where particles 
belonging to A^ distinct classes with hierarchical order diffuse as well interchange positions 
with rates depending on their relative hierarchy. We show that the exact solution can also 
be derived in the general case where the particles have arbitrary sizes. 
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Some extensions of our results can be made. A first and quite interesting generalization of 
our model happens when we allow molecules in any class to have size s = 0. Molecules of size 
zero do not occupy space on the lattice, having no hard-core exclusion effect. Consequently 
we may have, at a given lattice point, an arbitrary number of them. The Bethe-ansatz 
solution presented in the previous section is extended directly in this case (the equations are 
the same) and the eigenenergies are given by fixing in (|7l|)- ([73|) the appropriate sizes of the 
molecules. It is interesting to remark that particles of a given class c' (2, 3, ... , A^), with size 
Sc' = 0, contrary to the case Sc' > 1, where they "accelerate" the diffusion of the particles 
in classes c < c', now they "retard" the diffusive motion of these particles. The quantum 
Hamiltonian in the cases where the particles have size zero is obviously not given by (^ 
but can be written in terms of spin S = oo quantum chains. Another further extension of 
our model is obtained by considering an arbitrary mixture of molecules, where molecules 
in the same hierarchy may have distinct sizes. The results presented in [^ correspond to 
the particular case of this generalization where A^ = 1 (simple diffusion). For general A^ 



the S matrix we obtain in (^51) is also a solution of the Yang-Baxter equation (0), but 
the diagonalization of the transfer matrix of the associated inhomogeneous vertex model is 
more complicated. The Bethe-ansatz equations in the case of the asymmetric diffusion, with 



particles of unit size |1T0|JTT[| , or with arbitrary size |jT5[ , were used to obtain the finite-size 
corrections of the mass gap Gn of the associated quantum chain. The real part of these 
finite-size corrections are governed by the dynamical critical exponent z, i. e., 

Re(G'jv) ~ A^"'. (74) 

The calculation of the exponent z for the model presented in this paper, with particles of 
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arbitrary sizes, is presently in progress |30[ . 
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Figure Captions 



Figure 1 - Example of configurations of molecules with distinct sizes s in a lattice of size 
L = 6. The coordinates of the molecules are denoted by the black squares. 
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